arXiv:1508.02529vl [cond-mat.str-el] 11 Aug 2015 


Suppression of the impurity-induced local magnetism by the opening of a spin 

pseudogap in Ni-doped Sr 2 Cu0.3 

Yannic Utz, 1 Franziska Hammerath, 1,2 Satoshi Nishimoto, 1 Christian Hess, 1 Neela Sekhar Beesetty, 3 
Romuald Saint-Martin, 3 Alexandre Revcolevschi, 3 Bernd Buchner, 1,2 and Hans-Joachim Grate 1 

1 IFW-Dresden, Institute for Solid State Research, PF 270116, 01171 Dresden, Germany 
2 Institute for Solid State Physics, Dresden Technical University, TU-Dresden, 01062 Dresden, Germany 
3 SP2M-ICMMO UMR-CNRS 8182, Universite Paris-Sud, 91405 Orsay Cedex, France 

(Dated: August 12, 2015) 

The S' = 1/2 antiferromagnetic Heisenberg spin chain compound Sr 2 CuC >3 doped with 1% and 
2 % of Ni impurities has been studied by means of 63 Cu nuclear magnetic resonance. A strong 
decrease of the spin-lattice relaxation rate Tff 1 at low temperatures points toward a spin gap, 
while a stretching exponent A < 1 and a frequency dependence of Tff 1 indicate that this spin gap 
varies spatially and should rather be characterized as a spin pseudogap. The magnitude of the 
spin pseudogap scales with doping level. Our results therefore evidence the finite-size character of 
this phenomenon. Moreover, an unusual narrowing of the low temperature NMR lines reveals the 
suppression of the impurity-induced staggered paramagnetic response with increasing doping level. 

PACS numbers: 75.10.Pq, 75.40.Gb, 76.60.-k 


The one-dimensional (ID) S' = 1/2 antiferromagnetic 
Heisenberg model can be called the harmonic oscillator 
of quantum magnetism. Its integrability makes it exactly 
solvable and therefore it is often used as an archetype for 
low-dimensional quantum magnets in theory. Despite its 
simplicity, it shows a very unusual behavior. The ground 
state of this model is an example of an highly entangled 
many-body quantum state, which is characterized by a 
lack of long-range order even at absolute zero. Its ele¬ 
mentary excitations are exotic quasiparticle excitations 
with fractional quantum numbers, the 5 = 1/2 spinons, 
which can be excited with infinitely low energy, i. e. 
the excitation spectrum has no energy gap to the ground 
state [1-4]. Regardless of its fundamentally important 
role, it is difficult to find realizations of this model in 
nature. Small perturbations induced by impurities or 
by interchain interactions are expected to lead to gaps 
in the excitation spectra or to three-dimensional (3D) 
long-range ordering [5-10]. It is therefore important to 
perform clear-cut experiments which explicitly address 
these perturbations in a controlled manner. The inves¬ 
tigation of the staggered paramagnetic response around 
intrachain impurities has been proposed to be a valuable 
tool for this purpose [11]. 

The cuprate compound Sr 2 CuC >3 is known to be among 
the best realizations of the ID 5 = 1/2 antiferromagnetic 
Heisenberg model. There, the chains are realized by cor¬ 
ner sharing CuOq plaquettes with 5 = 1/2 on the cop¬ 
per site, which are mainly interacting along one crystal¬ 
lographic axis with a large exchange coupling of about 
J ~ 2000 K [12]. Weak static magnetism occurs only be¬ 
low T)v = 5.4 K [8, 13], which is low compared to the 
much larger exchange coupling J. However, recent stud¬ 
ies on doped variants of Sr 2 CuC >3 and the closely related 
double chain compound SrCu0 2 revealed the vulnerabil¬ 
ity of the originally gapless spinon excitation spectrum 


[14, 15] against the influence of impurities and disorder. 
63 Cu nuclear magnetic resonance and transport studies 
showed that doping Ca on the Sr site outside the chains 
breaks the integrability of the model and opens a spin 
gap of similar size in both compounds, which has been 
attributed to structural distortions and a concomitant 
bond disorder [16-19]. Inelastic neutron scattering dis¬ 
closed a striking impact of minor concentrations of intra¬ 
chain nickel impurities on the low-energy spin dynamics 
of the double chain compound SrCu0 2 [20]. The authors 
report the emergence of a spin pseudogap of the order of 
A « 90 K by replacing only 1 % of the 5 = 1/2 copper 
ions with 5 = 1 nickel impurities. Corresponding to their 
interpretation the nickel spin is fully screened. Therefore, 
the nickel ions effectively act as S=0 impurities and ba¬ 
sically cut the chains into segments with varying finite 
length l , which show finite-size spin gaps with magni¬ 
tudes proportional to l/l [6]. 

In this article, we show results on the doping- 
dependent effect of nickel impurities on the single chain 
compound Sr 2 Cu 03 based on nuclear magnetic resonance 
(NMR). Spin-lattice relaxation measurements reveal the 
opening of a spin pseudogap that scales linearly with the 
Ni content and therefore prove the finite-size nature of 
this phenomenon. However, the NMR spectra evidence a 
suppression of the impurity-induced staggered paramag¬ 
netic response with increasing impurity content, in stark 
contrast to what has been observed so far in gapped low 
dimensional spin systems [11, 21 -25]. 

The measurements were performed on high purity sin¬ 
gle crystals of Si^Cui-^Ni^Os (x = 0.01 and 0.02, la¬ 
beled Nil and Ni2 hereafter). The samples were prepared 
using the travelling solvent floating zone (TSFZ) method, 
since this compound undergoes a peritectic-type decom¬ 
position upon cooling from the melt [26]. The used start¬ 
ing powders of SrC0 3 , NiO and CuO were of 99.99 % pu- 
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rity. The crystals were grown at a growth rate of 1 mm/h 
under flowing oxygen atmosphere (50ml/min). A solvent 
pellet of composition 63 %CuO-37 %SrO was used to ini¬ 
tiate the crystal growth experiment. The crystals cleave 
readily along (h, 0,0). The high quality of the crystals 
has been checked by X-ray diffraction (phase determina¬ 
tion) and energy-dispersive X-ray spectroscopy (chemical 
composition) measurements. 

The 63 Cu NMR spectra have been measured with the 
standard Hahn spin echo method at a fixed frequency of 
80 MHz by sweeping the external magnetic field H and 
integrating the echo. At room temperature, the NMR 
spectra of the (/ = 3/2) 63 Cu nuclei consist of three 
narrow quadrupolar split lines with satellites which are 
only slightly affected by quadrupolar broadening [27] that 
results from the structural disorder mainly induced by 
the Ni-dopants. 

The 63 Cu NMR spin-lattice relaxation rate Tj -1 has 
been measured by the inversion recovery method on one 
of the quadrupolar split satellite lines for both doping 
levels. All T\ measurements were performed in magnetic 
fields close to 7T [28] with H parallel to the crystallo¬ 
graphic a axis. Accurate alignment of the samples has 
been achieved by utilizing the angle dependence of the 
second order quadrupolar shift of the 63 Cu main line. 
The recovery curves of the nuclear magnetization have 
been fit to the standard function for magnetic relax¬ 
ation of I = 3/2 nuclei measured on a satellite transition 
[29, 30]: 


M z {t) = M 0 1 


/(o.4e _(6t/Tl) * + 0.5e _(3t/Tl)> 
+ 0.le~^ T ^) 


(1) 


Mo is the equilibrium value of the nuclear magnetization, 
/ is ideally 2 for a complete inversion and the stretching 
exponent A < 1 accounts for a distribution of spin-lattice 
relaxation rates around a characteristic value T-/ 1 . For 
high temperatures, A ss 1 shows that there is a unique 
spin-lattice relaxation rate T) -1 . For lower temperatures, 
A < 1 indicates a spatial distribution of nuclei with dif¬ 
ferent spin-lattice relaxation rates. The ratio of T\ of 
the two Cu isotopes indicates purely magnetic relaxation 
over the whole temperature range. 

We will start our discussion with the resulting NMR 
spectra. Fig. 1 shows the 63 Cu main line for both dop¬ 
ing levels and various temperatures. Upon lowering the 
temperature, the spectra show a pronounced broadening. 
They develop shoulder structures and a splitting of the 
central peak. Both, the central line and the satellites, are 
equally affected (see appendices), which indicates that 
the broadening is of magnetic origin. This means that 
the spectra can be seen as a histogram of the distribu¬ 
tion of local magnetic fields. It is well known that such a 
magnetic broadening within antiferromagnetically corre¬ 
lated systems can be attributed to the presence of clouds 
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Figure 1: 63 Cu NMR main line of S^CuCU doped with 
1 % and 2 % of nickel at various temperatures obtained 
at a fixed frequency of 80 MHz by varying the external 
field. The intensity of the resonance lines is normalized 
to their maxima. The dashed lines indicate the 
expected I/Vt behavior of the shoulder feature and 
has been fit to its width at 100K. 



of field-induced staggered polarization around impurities 
[11, 21-25]. The shape and extension of such a local al¬ 
ternating magnetization (LAM) depends very much on 
the nature of the underlying spin system and the nature 
of its coupling to the impurity spin. 

A LAM has been observed previously in undoped 
S^CuCU, where it has been explained by open chain 
ends due to excess oxygen [31-33]. These chain ends 
break the translational invariance of the spin chain and 
lead to a local alternating susceptibility [XaitOz)]> which 
gives rise to a LAM in a magnetic field. It could be 
modeled based on the assumption of semi-infinite chains 
[34], which predicts a LAM with a maximum at a certain 
distance l = 0A8J/T [31] from the impurity and an ex¬ 
ponential decay for larger distances. Upon lowering the 
temperature, the maximum shifts further into the chain 
and increases with Xait,max oc 1/VT. In the NMR spec¬ 
tra, this causes a broad background with sharp edges, 
which broadens with decreasing temperature correspond¬ 
ing to AH oc 1 /y/T independent of the amount of chain 
breaks [31, 34], The intensity of the background should 
increase with decreasing temperature. 

We can identify the shoulder features as this broad 
background. Owing to a larger impurity concentration 
than in the undoped compound, the features are already 

























3 


well developed at 100 K for Nil. For Ni2, one can already 
observe shoulder features at 200 K. The dashed lines in 
Fig. 1 indicate the expected 1 /VT behavior [35]. The 
line shape clearly follows this trend. However, at lower 
temperatures the shoulder structure is smeared out. The 
onset of the smearing depends on the doping level. It sets 
in at higher temperatures for Ni2 than for Nil. Moreover, 
one can see that the resonance lines at low temperatures 
(below 100 K) are narrower for Ni2 than for Nil. This 
smearing of the shoulder structures and also the narrow¬ 
ing of the low-temperature resonance lines with increas¬ 
ing doping level are surprising and cannot be explained 
by the simple approach to Xait(aO mentioned above. 

At low enough temperatures, the LAM is predicted to 
extend over the whole chain segment and the assumption 
of the semi-infinite chain is no longer valid [9, 11, 33, 34], 
Due to the high impurity content, the finite size of the 
chain segments cannot be neglected in the investigated 
temperature range. However, it is not expected to result 
in a reduced width as compared to the 1/yT behavior or 
even a disappearance of the shoulder features. The devi¬ 
ations might be connected to the screening of the Ni spin. 
Since Ni is a magnetic impurity, an additional screening 
cloud is expected to contribute to the LAM [ 6 , 36]. More¬ 
over, susceptibility measurements on Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 
show that the nickel spin is in fact screened [37]. But 
such a screening cloud should rather enhance the width 
of the NMR line than suppress it. The suppression of the 
linewidth might also be connected to the results on the 
low-energy excitation spectrum, which will be discussed 
in the following. 

Fig. 2(a) shows the temperature dependence of T-f 1 
measured at the center of the resonance lines for Nil and 
M2. At high temperatures, Tf 1 is temperature inde¬ 
pendent for both dopings, as it is theoretically expected 
for antiferromagnetic S' = 1/2 Heisenberg chains [38, 39] 
and as it has been experimentally verified for the parent 
compound [15]. Below a certain crossover temperature, 
which is Tf n « 110 K for Nil and Tf i2 « 200 K for M2, 
Tf 1 shows a strong decrease by two orders of magni¬ 
tude toward low temperatures. The decrease of T-f 1 is 
accompanied by a decrease of the stretching exponent A 
[see Fig. 2(b)] and thus by a growing spatial distribution 
of spin-lattice relaxation rates, which levels off at lower 
temperatures. 

Due to the hyperfme coupling A± between nuclei and 
electrons, Tf 1 measures the imaginary part of the dy¬ 
namic spin susceptibility x" of the electronic spin system 
at the NMR frequency. For pure magnetic relaxation, it 
is given by 

( 2 ) 

— UJ 

Q 

On a more intuitive level, the relaxation mechanism 
can be described as the scattering of thermally excited 
spinons by the copper nuclei [40]. 
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Figure 2: 63 Cu NMR spin-lattice relaxation rate Tf 1 
(a) and stretching exponent A (b) of Sr 2 Cu 0 . 99 M 0 . 01 O 3 
(black circles) and Sr 2 Cu 0 . 98 Ni 0 . 02 O 3 (green squares) 
measured at the center of the high field satellites. Inset 
of (a): Arrhenius plot of the spin-lattice relaxation 
rates. 


Thus, the decrease in spin-lattice relaxation rates 
clearly indicates the depletion of low-lying states in the 
spin excitation spectrum, and therefore points toward a 
spin gap. However, the distribution of spin-lattice relax¬ 
ation rates, as indicated by A < 1, implies that this spin 
gap varies spatially and should rather be characterized 
as a spin pseudogap [41]. 

Usually, the magnitude of a spin gap is estimated by 
fitting the temperature dependence of Tf to an acti¬ 
vated behavior [16, 42-45] and using the activation en¬ 
ergy as an estimate for the spin gap. However, in our 
case, the spin-lattice relaxation rates do not decrease 
exponentially [see the inset of Fig. 2(a)]. This can be 
attributed to the spatial distribution of spin gaps, be¬ 
cause the fast relaxation stemming from nuclei exposed 
to small gaps will dominate the recovery process at low 
temperatures. We use the crossover temperature T* as 
an estimate for the average gap energy. T /) i2 is about 
twice as large as Tf n . Therefore, we conclude that the 
spin pseudogap is proportional to the doping level. This 
is in agreement with the assumption that the individual 
chain segments show gaps A a l/l and thus evidences 
the finite-size character of the spin pseudogap. The value 
of Tf n ss 110 K is close to the reported spin pseudo¬ 
gap A ss 90 K [20] of the double chain compound doped 
with 1 % of nickel, which suggests that the double chain 
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Figure 3: Spin-lattice relaxation rate Tf 1 (solid 
squares) and stretching exponent A (open circles) of 
Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 measured at different positions 
within the 63 Cu high field satellite (solid lines) at 
various temperatures in a magnetic field of 
fj>oH = 7.0493T with H\\a. The intensity of the 
resonance lines is scaled to arbitrary units. 


structure is not crucial to the gapping mechanism, simi¬ 
lar to what has been observed in the Ca-doped variants 
of SrCu 02 and Sr 2 Cu 03 [16, 18]. 

Having established the proportionality between the 
doping level and the average gap energy, we wanted to 
gain further knowledge about the spatial variation of spin 
gaps. Therefore, we investigated the frequency depen¬ 
dence of Tf 1 within the broad resonance lines. Fig. 3 
shows spin-lattice relaxation rates T-f 1 and stretching 
exponents A measured at different positions within the 
high field satellite of Nil [46]. The spectra themselves 
are also shown for guidance. In contrast to the spectra 
shown before, these spectra were obtained at a fixed field 
/j,qH = 7.0493 T by sweeping the frequency and adding 
the Fourier transforms of the echo signals (frequency 
step and sum method [47]). While Tf 1 is frequency- 
independent for T > 90 K, it shows a strong frequency de¬ 
pendence at lower temperatures. Spin-lattice relaxation 
rates at all positions decrease toward low temperatures, 
but the decrease of Tf 1 is less pronounced for larger dis¬ 
tances to the center of the resonance line. The frequency 
dependence of Tf 1 is accompanied by a frequency depen¬ 
dence of the stretching exponent A, which sets in already 


at T = 90 K. A is minimal at the center and larger at the 
outer parts of the resonance lines. This indicates that 
Cu nuclei which contribute to the outer parts of the res¬ 
onance lines probe a narrow distribution of small spin 
gaps, while Cu nuclei contributing to the center of the 
resonance lines probe a broad distribution of large and 
small spin gaps. As NMR is a spectroscopic method, we 
cannot distinguish if the gap differs only between chain 
segments of different length or if it also varies within in¬ 
dividual chain segments. If the shape of the LAM in 
real space would be known, such a distinction could be¬ 
come possible. Besides, we cannot exclude that addi¬ 
tional impurity-induced spin fluctuations the near chain 
ends lead to a variation of Tf 1 within single chain seg¬ 
ments and therefore contribute in combination with the 
LAM to the frequency dependence of T-f 1 . Such spin 
fluctuations may also enhance the deviation of Tf 1 from 
exponential behavior at low T [see the inset of Fig. 2(a)]. 
However, we can state that the largest gaps are measured 
at the center of the resonance lines and therefore by nu¬ 
clei not exposed to the LAM. This might be the key to 
understand the suppression of the LAM at low T and the 
frequency dependence of Tf . It suggests that the LAM 
is suppressed by the gap. This is supported by the fact 
that the gap increases with increasing impurity concen¬ 
tration, while the spectral broadening is reduced. In view 
of this considerable suppression of local magnetism, it is 
not surprising that the ordering temperature is strongly 
reduced too. We did not find any signature of magnetic 
ordering down to 4.2 K for both dopings. This is in agree¬ 
ment with recent measurements of the susceptibility and 
the specific heat of Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 , which do not show 
any transition down to 2K [37]. 

In summary, our 63 Cu NMR measurements on sin¬ 
gle crystals of the S = 1/2 spin chain S^Cui-^NfOs 
(x = 0 . 01 , 0 . 02 ) show a strong impact of minor concentra¬ 
tions of nickel on the low-energy spin dynamics and the 
local susceptibility of the spin chains. We find a doping- 
dependent spin pseudogap behavior, which evidences the 
finite-size character of this phenomenon. The NMR spec¬ 
tra show a local alternating magnetization around the 
nickel impurities. Its suppression at low temperatures, 
the variation of Tf within the broad resonance lines 
and the strong reduction of the ordering temperature 
are most probably consequences of the spin pseudogap, 
which reduces low-energy antiferromagnetic fluctuations. 
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Appendices 

Attempt to model the distribution of Tj -1 to fit the 
recovery curves 

As demonstrated in the main article, the measured re¬ 
covery curves of the nuclear magnetization at low temper¬ 
atures are stretched due to a distribution of spin lattice 
relaxation rates Tf 1 . We took this distribution into ac¬ 
count by fitting the phenomenological stretched exponen¬ 
tial function to the measured recovery curves. However, 
this approach just gives a rough idea about the width 
of the distribution in terms of the stretching exponent 
A. It is also possible to convolute the recovery function 
M z (t, T\) with a chosen probability distribution function 
of T\ . In accordance with the approach described by [20], 
we modeled the distribution of T\ as result of the distri¬ 
bution of gaps due to a given distribution of chain length. 
Following [20], the probability to find a copper atom in¬ 
side a non-interrupted segment of length L for a totally 
random distribution of defects with a concentration x is 
given by 

p. _ L-7?{l-x) L _ 2 xL-r 

p TZi L - x 2 i l ~ x ) L ~ L{ ' 3 

A segment of length L has an energy gap A L = A 0 /L 
with Ao = 3.65J [ 6 , 20]. Assuming an activated behav¬ 
ior, the temperature dependence of T\ for copper nuclei 
within a chain segment of length L is given by 

Ti(T,L) = Ti i0O exp(A L /T) = T 1>00 exp(A 0 /T/T). (4) 

Tj oc labels the high temperature limit, which is inde¬ 
pendent of thegap magnitude. The recovery function for 
magnetic relaxation of / = 3/2 nuclei measured on a 
satellite transition with a unique spin-lattice relaxation 
time Tj is [29]: 


M. 


: (t, Tj) = M 0 1 - /(o.4e- w/Tl + 0.5e _3 * /Tl (5) 

+ 0 . 1 e" t/Tl ") 


The recovery function based on the distribution of 
chain length can then be obtained as weighted sum over 
L of Eq. (5), where Tj is given by Eq. (4): 


M z ( t ) = Y, M z (*, Ti (T, L)) ■ P L . ( 6 ) 

l 

For the fitting, we set A 0 = 3.65 • 2000 K = 7300 K and 
T corresponding to the current temperature fixed, x, 
Tqoo, Mo, and / were used as adjustable parameters. 

However, it was not possible to obtain satisfactorily 
fits. Fig. 4 shows a curve measured at 25 K and its fit 
as an example. One can see, that the fitting curve does 



Figure 4: Recovery curve of the nuclear magnetization 
measured at the center of the 63 Cu high field satellite of 
Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 at 25 K in a magnetic field 
HoH = 7.0493 T oriented parallel to the crystallographic 
a-axis. The green line illustrates the fit with the 
described model. The fitting parameters are: 
x = 0.020 ± 0.003, Ti )00 = (0.22 ± 0.07) ms, 

M 0 = 1.36 ± 0.08 and / = 1.66 ± 0.06. T = 25 K and 
A 0 = 7300 K where set fixed during the fitting 
procedure. The blue line is simulated corresponding to 
the model with the same parameters as the fit except of 
x = 0.01. The red line depicts a stretched fit for 
comparison, with Ti = (7.5 ± 0.3) ms and 


A = 0.47 ±0.01. 


not reach full recovery for very long time in contrast to 
the measured curve. Moreover, the resulting doping level 
x = 0.02 is double as large as the nominal one. Moreover, 
the parameters x and Ti ;(X) vary a lot over temperature, 
which should not be the case. Setting x and Ti >00 to 
reasonable values and holding them constant during the 
fitting procedure resulted in even worse fit quality. 


The reason for the failed fitting could be that the orig¬ 
inal distribution of gaps is spread over the broad reso¬ 
nance lines, so that only excerpts of the original distri¬ 
bution are measured at different positions, as stated in 
the main article. This is supported by the differences be¬ 
tween measured and fitted curves. In Fig. 4 one can see 
that the measured curve realizes a narrower distribution 
of T\ than assumed by the model. This is manifested by 
the fact that the initial recovery of the measured curve is 
much slower than it should be for x = 0.01 (see simulated 
curve in Fig. 4) and that the very slow recovery at long 
times it not visible in the measured data. 
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Spectra and Linewidth 

Fig. 5 and 6 show field-swept Cu NMR spectra of 
Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 and Sr 2 C 1 io. 98 Nio. 2 O 3 obtained with 
the standard Hahn spin echo method at a fixed frequency 
of 80 MHz by sweeping the magnetic field H along the 
crystallographic a axis and integrating the echo. They 
contain all three resonance lines of the 63 Cu spectrum 
and the 65 Cu high field satellite. 

As the shape of the resonance lines changes drastically 


over temperature, we use the second moment as a mea¬ 
sure of linewidth (see Fig. 7 for the linewidth of the 63 Cu 
main line). For temperatures T < 25 K, the tails of the 
main line have an overlap with the tails of the neigh¬ 
boring satellite lines. In these cases, we integrated only 
between the two minima surrounding the mainline. This 
leads to incorrect values for temperatures T < 25 K. In 
spite of these limitations, Fig. 7 shows that the increase 
of the linewidth towards low temperatures is smaller for 
2 % Ni doping than for 1 % Ni doping. 



Echo Intensity (arb. u.) Echo Intensity (arb. u) 
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Figure 6 : Field-swept Cu NMR spectra of Sr 2 Cuo. 9 sNio. 2 O 3 . 
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Figure 7: Linewidth a of the 63 Cu main line within the spectra shown 
above, determined as the square root of the second moment, as a function 
of temperature. The solid lines are guides for the eyes. The dashed line 
indicates the temperature below which the second moment is not a reliable 
measure of the line width any more due to the overlap with the satellite 

lines. 











